a-Scale Decoupling of the Mechanical Relaxation and Diverging Shear Wave 
Propagation Lengthscale in Triphenylphosphite 
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We have performed depolarized Impulsive Stimulated Scattering experiments to observe shear 
acoustic phonons in supercooled triphenylphosphite (TPP) from ~10 - 500 MHz. These measure- 
ments, in tandem with previously performed longitudinal and shear measurements, permit further 
analyses of the relaxation dynamics of TPP within the framework of the mode coupling theory 
(MCT). Our results provide evidence of a coupling between the shear and longitudinal degrees of 
freedom up to a decoupling temperature Tc — 231 K. A lower bound length scale of shear wave 
propagation in liquids verified the exponent predicted by theory in the vicinity of the decoupling 
temperature. 

PACS numbers: 64.70.ph, 61.20.Lc, 43. 35. Mr, 81.05.Kf 



I. INTRODUCTION 

The most prominent characteristic of the frequency re- 
laxation spectrum of glass forming liquids is the a peak, 
a broad, yet distinct feature whose origin lies in the col- 
lective degrees of freedom which undergo kinetic arrest 
during vitrification [ll-Q. As the transition from liquid 
to glass occurs either upon cooling or the application of 
high pressure, the characteristic time scale of a relax- 
ation varies by 16 orders of magnitude, from the ■~100 fs 
timescale associated with an attempt frequency in the 
liquid state [4] to the 100 s timescale (arbitrarily) cho- 
sen at Tg (3|]. One of the fundamental tasks presented 
to the spectroscopy of glass forming liquids is thus the 
complete characterization of the a relaxation across this 
span of frequencies. 

To date, this challenge has best been met by dielec- 
tric spectroscopy, whose wide dynamical range remains 
unparalleled in the study of glass forming liquids [B-Q]- 
Nevertheless, dielectric spectroscopy of molecular liquids 
probes mainly orientational relaxation, while the glass 
transition is more naturally understood in terms of the 
degrees of freedom associated with fiow, i.e. the den- 
sity, which may be probed by longitudinal acoustic waves, 
and the transverse current, which couples to shear waves 

M- 

Despite the importance of both of these degrees of free- 
dom, the bulk of the literature on structural relaxation 
has focussed on study of the longitudinal modulus due to 
the relatively limited means by which shear waves may 
be generated and probed. Acoustic transducer [l^, ITl| . 
Impulsive Stimulated Scattering (ISS ) fl2| - [l3 | and high- 
frequency forced Brillouin Scattering [15| techniques can 
continuously cover the longitudinal spectrum up to the 
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mid- GHz regime. However, broadband shear wave gener- 
ation is more accessible experimentally at lower frequen- 
cies and exhibits a gap in the range from the ~ 10 MHz 
limit of transducer methods (lol . [ll| to the low GHz fre- 
quencies probed by depolarized (VH) Brillouin Scatter- 
ing |16l - l21| . The vast majority of depolarized Brillouin 
Scattering studies are conducted in a single scattering 
geometry, permitting observation of only one acoustic 
wavevector. Thus, relaxation information must be ex- 
tracted by complex modeling incorporating contributions 
from the various channels responsible for the depolarized 
scattering of light. Obtaining spectral information from 
a collection of acoustic frequencies would allow a direct 
determination of the shear relaxation spectrum without 
the need for multiparameter fitting. 

Recent advancements in the generation of narrowband 
shear waves [2^ using the technique of Impulsive Stim- 
ulated Brillouin Scattering (ISBS) in the depolarized 
geometry [2^ now permit study of the shear a relax- 
ation spectrum across a broad range of frequencies. In 
this study, we employ ISBS to focus on tests of a first- 
principles theory of t he g lass transition, the mode cou- 
pling theory (MCT) |24l - [26| . While the natural vari- 
ables of the theory are the density fluctuations that con- 
stitute the longitudinal modes, MCT suggests relation- 
ships among the temperature dependencies of the var- 
ious relaxing variables (e.g., transverse current fluctu- 
ations, density fluctuations, orientational fluctuations, 
etc.) through the property of a-scale coupling [26|. This 
prediction states that if two variables, A and B, couple 
to density fluctuations, their associated relaxation times, 
ta{T) = Cat{T) and tb{T) = Cbt{T), are equivalent, 
up to a temperature independent factor, to a universal 
relaxation time t{T) characteristic of the main a relax- 
ation process. Decoupling may then occur at the MCT 
critical temperature Tc when the relationship between 
these relaxation times can break down. Prior experimen- 
tal studies that have examined a-coupling have focused 
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on, e.g., the coupling of dielectric to rheological variables 
[27| , dielectric to shear mechanical relaxation [28| , and of 
translational to rotational diffusion [2I). However, there 
have been no studies that directly link elastic degrees of 
freedom with each other. 



MCT predictions have also been formulated concern- 
ing the transverse current fluctuations in their own right, 
in particular the power-law divergence of an upper bound 
length scale for shear wave propagation derived for hard 
spheres in a Percus-Yevick approximation [s^. While 
such upper bound has been established in the literature 
Q in terms of a wavelength at which the correspond- 
ing shear acoustic frequency goes to zero, this particular 
scaling relationship with temperature is unique to MCT 
and, to our knowledge, has not undergone scrutiny in the 
lab. 



Below, we provide a test of these assertions by building 
upon prior work [3l|, , and present an expanded study 
of the shear acoustic behavior of the fragile glass-former 
triphenylphosphite (TPP) (Tg = 202 K ;.33]). Taken to- 
gether, these experiments fill in the gap between 10 MHz 
and 1 GHz in the shear relaxaiton spectrum, and thus 
constitute the broadest bandwidth shear acoustic mea- 
surements performed optically any glass forming system 
to our knowledge. Our measurement approach thus en- 
ables characterization of shear relaxation in glass-forming 
liquids in this frequency regime. Although supercooled 
liquid a relaxation spectra are typically broader than two 
decades at any temperature and, as T is varied, the spec- 
tra move across many decades, the 10-1000 MHz range 
is sufficient to permit assessment of whether shear relax- 
ation dynamics are consistent with relaxation dynamics 
determined from longitudinal wave measurements cover- 
ing a comparable spectral width at common tempera- 
tures. They are also sufficient to test predictions about 
the low-frequency limit of shear wave propagation as a 
function of temperature. 



We begin with a general overview of depolarized ISBS. 
This motivates a derivation of the signal generated in an 
ISBS experiment geared toward measuring shear waves 
from simplified versions of the hydrodynamic equations 
of motion. After a brief review of the experimental appa- 
ratus used, we present data collected in the temperature 
range from 220 K to 250 K. Our measured shear acoustic 
frequencies show significant acoustic dispersion, indica- 
tive of the complex relaxation dynamics characteristic 
of supercooled liquids thus providing a means of test- 
ing the MCT predictions described above. In particular, 
we see evidence of a-decoupling, yielding an estimate of 
the MCT Tc of 231 K. Using this value of the crossover 
temperature yields a power law for the lower limit of 
shear wave propagation in line with predictions derived 
in Ref . [30] . These results are discussed within the frame- 
work of the MCT. 



detector 




FIG. 1. Schematic illustration of the ISS setup. Both the 
pump and probe beams are incident on the PM and their ±1 
diffraction orders are recombined at the sample. In the case of 
a depolarized experiment, half-waveplates in the path of each 
of the four beams are used to create a polarization grating 
pattern in the sample plane, as shown in the figure inset. 



II. ISBS - GENERAL CONSIDERATIONS 

In a typical ISS experiment, conducted in a hetero- 
dyned four-wave mixing geometry, light from a pulsed 
laser is incident on a diffractive optical element, typically 
a binary phase mask (PM) pattern, and split into two 
parts (±1 diffraction orders; other orders are blocked) 
that are recombined at an angle 9 as depicted in Fig. [T] 
The crossed excitation pulses generate an acoustic wave 
whose wavelength A is given by 



A 



2 sin 61/2 



(1) 



where A is the excitation laser wavelength. Probe light 
(in the present case from a CW diode laser) is also inci- 
dent on a phase mask pattern (the same one or another 
with the same spatial period) and split into two parts 
that are recombined at the sample to serve as probe and 
reference beams. The signal arises from diffraction of 
probe light off the acoustic wave and any other spatially 
periodic responses induced by the excitation pulses. The 
diffracted signal field is superposed with the reference 
field for heterodyned time-resolved detection of the sig- 
nal, which typically shows damped acoustic oscillations 
from which the acoustic frequency and damping rate can 
be determined. 

The polarizations of the beams in the present measure- 
ments are all vertical (V) or horizontal (H) relative to the 
scattering plane. We denote the polarizations of the exci- 
tation fields, probe field, and reference/signal field in that 
order, i.e. VHVH denotes V and H excitation, V probe, 
and H reference and signal polarizations while WW de- 
notes all V polarizations. All the measurements reported 
herein were conducted in the VHVH configuration. 

In a depolarized (VHVH) experiment, each of the 
pump arms carries a different polarization and the en- 
suing grating is described as an alternating polarization 
pattern, depicted in the Fig. [T] inset. It is the regions 
of linear polarization that perform electrostrictive work, 
deforming the excited region in a fashion that generates 
counterpropagating shear acoustic waves with a driving 
force that scales with wavevector magnitude q = 27r/A 
[23} . The driving force also scales with a geometrical 
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factor of cos 0/2 due to a diminishing of the contrast ra- 
tio of the polarization grating with increased angle. In 
the small 6 limit, the signal from shear waves becomes 
markedly stronger as the wavevector is increased. In this 
case, the diffracted signal polarization is rotated 90° from 
that of the incident probe light, analogous to depolarized 
(VH) Brillouin scattering. Finally, we note that the ex- 
citation pulses may also induce molecular orientational 
responses that can contribute to signal, analogous to de- 
polarized quasielastic scattering [3l|, I34l - l39l | . 



III. ISBS - THEORETICAL BACKGROUND 

A full treatment of signal in an ISS experiment includes 
longitudinal and transverse acoustic modes, molecular 
orientations, and other nonoscillatory degrees of freedom 
IssUioj . Here, we present a slightly simplified deriva- 
tion that allows us to understand the presence of shear 
and orientational contributions to the signals. 

We start with the generalized, linearized hydrody- 
namic equations of motion 41]. First is the momentum 
density conservation law. 



^pu(r,t) 



V • n(r,i) = 



(2) 



where p is the density, u(r, t) the velocity, and n(r, t) the 
stress tensor. The expression for the stress tensor reads 



+'5a,3(3?7-C)V- u (r,t) 

p,Qa[3{Y,t) + Fai3{T,t). 



(3) 



Here, S^p is the Kronecker delta, P{r,t) is the pressure, 
rj is the shear viscosity, ^ is the longitudinal viscosity, 
Qa/3(r,t) is the orientational variable, fi expresses the 
coupling of translational force due to rotational motion, 
and r is the spatial coordinate. Fnfiiy, t) is the external 
shearing stress of the laser field [ifl, E^l , assumed to be 
temporally impulsive and spatially periodic, i.e., Fap — 
S{t) cos q • r, and Tap is the rate of strain, defined by 



dua{r,t) dup{r,t) 



drt. 



(4) 



The orientational variable obeys its own equation of 
motion, given by 



dQap[v,t) 
dt 



-TrQ^p {r,t)+ £,Tap {Y,t) + Qo{v,t) (5) 



where F^ is the orientational relaxation rate, ^ is the 
torque due to translational motion, and (3o(r,i) is the 
torque exerted by the laser [20|, HI] . As with the laser- 
induced stress, the torque will also be modeled as tem- 
porally impulsive and spatially periodic. This equation 
of motion assumes Debye orientational relaxation as a 
computational convenience. 



We set the grating wavevector in the x direction and 
the transverse direction to be y. Since we are interested 
in shear waves, we select the transverse elements of the 
above equations, and after a Fourier-Laplace transform 
defined as 



-^^{/(r,t)} 



we are left with 



dt 



drf(r,t)e''i' 



(6) 



psuy{q,s) = -iqn^^iq, s) 
U^y{q,s) = -rjsT^y{q,s) + 
sQxyici, s) = -TRQxy{q, s) 



(7) 

tJ-Qxy{(l,s) (8) 

^r-y{q,s)+Qo. (9) 



T^xy 
^0 



For simplicity, we assume that the measured signal is 
proportional to the molecular polarizability anisotropy 
[42| , i.e. to Qap , so we solve for the orientational variable 
to yield 



Qoisp + q^T]) 



{s + TR){3p + q^r]) -Cw^ 



(10) 



In order to arrive at an analytic solution, we make the ap- 
proximation that ^p.q'^ is a coupling of higher order that 
can be ignored in the solution of the equations of mo- 
tion. Physically, we are making the approximation that 
we may ignore the recoupling of the orientational degree 
of freedom to itself via rotational-translational coupling. 
In this approximation, the expression for Qxy separates 
as 



Qxyiq^s) = 



xy 



(s + ^r) [sp + q^v) •'^ + T/ 



(11) 



In order to extract the effect of structural relaxation 
dynamics on signal, we model the a peak by Debye relax- 
ation as 



, where Coo is the infinite frequency speed of 
sound and Tg is the characteristic shear relaxation time. 
This, in effect, makes rj complex, which enables an elas- 
tic component of the shear response to emerge from the 
equation of motion[3] 41]. In this model, the above equa- 
tions can be solved for Q^y to yield 



2 T^xy 



1 



Qo 



(12) 



Equation (fT2|) can be recast into the form: 

^q'F^y 2r, + s 



Qxy(q, s) 



P (s + Ffl)(s + r, + iw,)(s + r, -zw,) 



s + Tr 

where the shear acoustic damping rate Fs is given by 



(13) 



1 
2^ 



(14) 
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and the frequency of oscillation lOs by 




2r, 



(15) 



The acoustic damping rate T g thus comprises the struc- 
tural relaxation dynamics. We also note that Us may go 
to zero for finite wavevector when 



1 



9o 



(16) 



Separation of equation [T3] by partial fractions yields a 
time-domain solution 

Qxy{q,t) = ii— 2-[Ae-i^»*sin(u;st) 

+B (e"^«* - e-^=*cos(w,t))] 



where 



A 



and 



B 



u1 + r^r, - 

{Tg-Tnf+uji) 



(2r, 



((r.-r«)' 



(18) 



(19) 



The solution represented by equation [IT] comprises two 
pieces. The term proportional to q^F^^ is due to the 
shear acoustic response, and the other term proportional 
Qo is a decaying exponential independent of q. This ori- 
entational response is the optical kerr effect (OKE) sig- 
nal. We also note the presence of an orientational con- 
tribution to signal in the acoustic response, due to the 
rotational-translational coupling. 

In order to derive the relaxation spectrum from 
our data, we consider a frequency-dependent modulus 
G*{s] = G'{s) + iG"{s) which obeys the dispersion rela- 
tion (4l 



ps^ + G{s)q'^ = 0, 



(20) 



and yields the following expressions for the real and imag- 
inary parts of the shear modulus from the acoustic fre- 
quency and damping rate 



q^ 



G [us) = p- 



q^ 



(21) 
(22) 



As equation[5n]has been derived considering the strain, in 
order to compare it with the results of the above analysis, 
we must solve for the strain from the original equations 
of motion. This gives the dispersion relation 



Comparison between equations [20l and [23l vields the con- 
nection between the elastic modulus and the relaxation 
spectrum r/{q,LUs) 



G'iujs) 



G"{ujs 



= -uJslm[r]{q,uJs)] 
= ujsRe[ri{q,ujs)]. 



(24) 
(25) 



ps 



■qsq = 0. 



(23) 



We conclude by noting that orientational responses of 
anisotropic molecules can be induced not only by the ex- 
citation pulses, as in the case of OKE, but also by flow 
that occurs due to the induced density changes [i^l ■ Both 
of these sources lead to signals that can be suppressed by 
proper selection of probe and signal polarizations. This 
step was impractical when measuring shear responses, as 
a choice of polarization which would suppress the orien- 
tational response often reduced the already weak shear 
signal beyond the limit of detection. Eliminating this 
contribution was also deemed unnecessary, as the aim 
of this study was to directly probe the shear relaxation 
spectrum via narrowband acoustic measurements of fre- 
quencies and damping rates. 



IV. ISBS ~ EXPERIMENTAL DETAILS 

The pump laser system used for these studies was an 
Yb:KWG High-Q FemtoRegen lasing at 1035 nm and 
producing pulses of 500 /iJ at a repetition rate of 1 kHz, 
although 150 /iJ was routinely used to avoid cumulative 
degradation of the sample. While a 300 fs compressed 
pulse width FWHM was typical, we bypassed the com- 
pressor to retrieve pulses directly from the regen that 
were 80 ps in duration in order to avoid sample damage 
at high peak powers, yet remain in the impulsive limit 
relative to the oscillation period. The excitations beams 
were cylindrically focussed to a spot that was 2.5 mm in 
the grating wavevector dimension and 100 puv in the per- 
pendicular dimension so that the acoustic waves would 
have many periods and the decay of signal would be due 
primarily to acoustic damping rather than propagation 
away from the excitation and probing region of the sam- 
ple. 

As a probe, we used a Sanyo DL8032-001 CW diode 
laser output at 830 nm with 150 mW power focused to a 
spot of 1 mm in the grating dimension by 50 /im in the 
perpendicular dimension. We also used a single phase 
mask optimized for diffraction into ±1 order at 800 nm 
for both pump and probe beams, and we utilized two-lens 
2:1 imaging with Thorlabs' NIR achromats to recombine 
the beams at the sample. The local oscillator was at- 
tenuated by a factor of lO^'^. Approximately 30% of the 
pump power was lost into zero order with this configu- 
ration, but the pump intensity still had to be reduced 
significantly to avoid unwanted nonlinear effects. 

In order to generate the polarization grating, we in- 
serted A/2 waveplates into each of the beams. These 
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waveplates were held in precision rotation mounts to pro- 
vide accurate alignment of the relative polarizations of 
the V and H polarized beams. This set the upper limit 
on the grating spacing that could be achieved in our mea- 
surements - for longer wavelengths, the beams came close 
enough together to be clipped by the rotation mounts. 
This issue was addressed by imaging with longer focal 
length optics. The signal was collected in a Cummings 
Electronics Laboratories model 3031-0003 detector and 
recorded by a Tektronix Model TDS-7404 oscilloscope. 
The shear signals were weak and required 10,000 aver- 
ages, resulting in data acquisition times of a few minutes. 

TPP at 97% nominal purity was purchased from Alfa 
Aesar and had both water and volatile impurities re- 
moved by heating under vacuum with the drying agent 
MgS04 immersed in the liquid. The sample was then 
transferred to a cell with movable windows [31 via fil- 
tering through a millipore 0.22 /im teflon filter. Af- 
ter loading, the cell was placed in a Janis ST-IOO-H 
cryostat where the temperature was measured with a 
Lakeshore model PT-102 platinum resistor immersed di- 
rectly within the liquid, and monitored and controlled 
with a Lakeshore 331 temperature controller. 

The grating spacings examined in this study were 
2.33 /xm, 3.65 /im, 6.70 /xm, 7.61 /Ltm, 9.14 /im, 10.2 yum, 
11.7 /xm, 13.7 /xm, 15.7 fim, 18.3 ^m, 21.3 /xm, 24.9 /xm, 
28.5 /Ltm, 33.0 /xm, 38.1 /xm, 42.6 /xm, and 50.7 /Ltm, while 
data for 0.48 /xm, 1.52 /xm, 3.14 /xm, and 4.55 /xm grating 
spacings were taken from prior reported results [sil, [s^] . 
The acoustic wavelength was calibrated through ISS mea- 
surements in ethylene glycol, for which the speed of sound 
is known to a high degree of accuracy [43l|. When the 
samples were cooled to the desired temperature, the cool- 
ing rate never exceeded 6 K/min, with 2 K/min being 
typical. Data were taken at fixed wavevector every 2 K 
from 220 K to 250 K upon warming as we found crys- 
tallization was less likely to occur upon warming than 
cooling. Only a few measurements could be taken with- 
out having to thermally cycle the liquid, as it invariably 
crystallized. We noticed that the tendency toward crys- 
tallization was particularly pronounced in the tempera- 
ture range between 234 K and 242 K. After a few days 
of use, the sample was observed to develop a slightly 
cloudy yellowish hue, and so was replaced by a new one. 
The yellowish samples tended more readily toward pho- 
toinduced damage, as well as crystallization, than the 
original, clear samples. Comparison of the signals ob- 
tained from the degraded samples and fresh ones yielded 
the same frequency and damping rate values, indicating 
that uncertainties in either of these quantities were due 
mainly to noise in the data. 



V. RESULTS AND DISCUSSION 

The results of several VHVH experiments performed at 
A=2.33 /im grating spacing are shown in Fig. [5J There 
is an initial spike due to the non-resonant electronic re- 
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FIG. 2. Shear waves in TPP with A = 2.33 /xm at (a) 220 K, 
(b) 234 K, and (c) 246 K. The data are in black and the fits 
in red. As the temperature is increased, the acoustic wave 
becomes more heavily damped. At higher temperatures, we 
also note the presence of the orientational relaxation, which 
appears to skew the signal such that the oscillations do not 
occur about the zero baseline. 



spouse. Immediately following this hyperpolarizability 
peak are oscillations about the zero baseline from the 
counterpropagating shear waves. At a sample tempera- 
ture of 220 K, these oscillations are seen to disappear on 
the scale of tens of nanoseconds due to acoustic damping. 
As the sample is warmed, the frequency is observed to 
decrease and the damping to increase dramatically. At 
sufficiently high temperatures, the shear wave becomes 
overdamped. We also note that at some temperatures, 
the signature of orientational relaxation is observed as a 
nonoscillatory decay component in the signal. 

Another illustration of the influence of relaxation dy- 
namics on the signal may be obtained by examining data 
from a collection of wavevectors at a common tempera- 
ture, as depicted in Fig. [3] where we provide data recorded 
with 10.2 /im, 21.3 /xm, and 44.2 /tm grating spacings 
at 220 K. Data with a fourth wavelength, 2.33 /tm, are 
shown in Fig. [2^. As the wavevector and the frequency 
are reduced, the acoustic oscillation period increases to- 
ward the characteristic relaxation timescale and there- 
fore the shear wave is more heavily damped. The signals 
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FIG. 3. Shear acoustic signal in TPP at 220 K for (a) 10.2 ^m, 
(b) 21.3 pim, and (c) 44.2 /xm grating spacings. The data are 
in black, and the fits are in red. We note that already by 
220 K at 44.2 ^m, only a few acoustic cycles are observed, 
and that the signal intensity is reduced. 



at larger grating spacings are weaker due to the linear 
q-dependence of the excitation efficiency. 

Based on the analysis of Sec. IIIIl time-domain signals 
were fit to the function 

I(t) = A' exp(-r5t) 'sm(ujst + (p)+B' expi-Vnt) + Cd{t) 

(26) 

which was convolved with the instrument response func- 
tion, modelled here by a Gaussian with duration 0.262 ns. 
The convolution was necessary to determine the true 
t = for the experiment. Here, A' is the acoustic ampli- 
tude, Ts is the shear damping rate and los is the shear fre- 
quency. (/) is a phase which accounts for the cosine term in 
Eq. 1171 which only becomes important when the damping 
is strong. In the next term, B' is the optical Kerr effect 
signal amplitude and is the orientational relaxation 
rate, and in the last term C is the strength of the hyper- 
polarizability spike. As discussed in Sec. Illli it is a sim- 
plification to model the orientational behavior by a sin- 
gle decaying exponential js^]. A more accurate descrip- 
tion might be in terms of a Kohlrausch- Williams- Watts 
stretched exponential function e(~*/*'"'») (0 < /? < 1), as 
is used commonly for time-domain relaxation; however. 
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FIG. 4. Shear speed of sound in TPP as a function of tempera- 
ture for a variety of grating spacings. We note that the highest 
temperature for which we could observe shear waves increases 
with the decrease of grating spacing. Data at A = 1.52 jj,m 
are from l3l|,[3|]. 



since the orientational signal contributions are weak, we 
were able to obtain excellent fits with fewer parameters 
using a single exponential form. 

The obtained values of the the shear acoustic veloc- 
ity Cs = ujs/q at a collection of wavevectors are shown 
in Fig. m while in Fig. [S] the scaled damping rates are 
shown. Both figures incorporate the data from Refs. ^31} 
and (32I . Data at other wavevectors were consistent with 
those shown, but are omitted from this and further plots 
for clarity. Two features of the data are immediately 
evident in these plots: first, we observe significant acous- 
tic dispersion for the shear waves across all temperatures 
studied, and this dispersion increases dramatically as the 
temperature is raised. The second feature we note is 
that, as a result of the dispersion and the shear soften- 
ing it represents, at each temperature above 240 K there 
is a wavelength above which we are unable to observe 
the shear wave in our measurements due to its increased 
damping and reduced signal strength. This wavelength 
is observed to decrease as temperature increases. 

From the fitted values for the shear frequency lOs and 
the damping rate F^, we may compute a value of the 
reactive and dissipative shear moduli from equations [21] 
and mi The density values we obtained from data in 
reference (45| were fit to a quadratic function as 



p{T) ^ 1.507 [g/cm^ 



-1.3x10"^ T [Kl 



-6. 



ixlO^^T^ 



(27) 

Figures H] and [7] show plots of the real and imaginary 
parts of the shear modulus, respectively, as functions of 
temperature. These plots are for the same collection of 
wavevectors for which we have plotted the velocity and 
damping information. As in the plot of the velocity, we 
see the softening of the modulus at higher temperatures. 
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FIG. 5. Scaled shear damping rate in TPP as a function of 
temperature for a variety of grating spacings. Longer wave- 
length acoustic waves have a higher scaled damping rate which 
increases with temperature. Data at A = 1.52 /^m are from 

MM- 



The imaginary part shows generally monotonic behav- 
ior as a function of temperature as well, except for the 
1.52 /im and 10.2 fxm data, which show a small decrease 
in the imaginary part at higher temperatures, a feature 
which is only weakly evident in the damping rate itself 
represented in Fig. [S] This is likely due to the oscillation 
period of our shear wave exceeding the characteristic re- 
laxation time Ts , permitting observation of a piece of the 
low-frequency side of the relaxation curve. We generally 
did not observe this trend at most wavevectors, as the 
shear wave signal became either too weak or too strongly 
damped to be observed. 

The moduli at each temperature were plotted as a func- 
tion of frequency and then fit to the Havriliak-Negami 
relaxation function 



G*K) = Go 



1 - 



1 



(1 + {iiO.TsYf 



(28) 



in order to extract the shear relaxation spectrum (Go 
for all temperatures, by definition). Since time- 
temperature superposition has been observed to hold for 
a relaxation in many liquids and in TPP (at lower tem- 
peratures than measured here) [ioj . all spectra were fit 
simultaneously with the spectral parameters a and /3 act- 
ing as temperature-independent global variables and r^. 
The value for the temperature dependent infinite shear 
modulus Goo = pc^ was taken from depolarized Bril- 
louin scattering measurements performed by Chappell 
and Kivelson [47[ in a linear extrapolation to colder tem- 
peratures as 

Coo = 2620.4 [m/s] - 7.93T [K] (29) 
where the above expression has been obtained by using 
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FIG. 6. G'(r) at a number of grating spacings. As the tem- 
perature is increased, the shear modulus is observed to soften 
considerably. Again, the highest temperature for which we 
could observe shear waves increases with the decrease of grat- 
ing spacing. Data at A = 1.52 /im are from [Sll. |3^. 
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FIG. 7. G"{T) at a number of grating spacings. For most 
wavevectors examined in this study, we only observed a rise 
in the value of G"{T) with temperature. At a handful of 
wavevectors, we were able to see the low-frequency side of the 
relaxation peak through, as is visible for 27r/o = 1.52 ^m and 
10.2 fim. Data at A = 1.52 /im are from [3ll. l32|. 



polarized (VV) Brillouin scattering data from longitudi- 
nal acoustic phonons to derive a temperature dependent 
refractive index ^4^] in combination with the shear data 
of the reference. We note that these data fall on a com- 
mon line with our lower-T shear data of Fig. [6] at the 
highest frequencies, i.e., the shortest wavelengths (A = 
1.52 /im and 2.33 /im). 

Three representative plots of the complex shear modu- 
lus with corresponding fits are shown in Figs.[S^ -[5}:. As 
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FIG. 8. Plots of the real (G') and imaginary (G") moduli 
of TPP at (a) 224 K, (b) 234 K, and (c) 242 K. At each 
temperature, a different segment of the shear acoustic relax- 
ation spectrum is present within the experimental frequency 
window. 



the temperature is increased, the shear relaxation spec- 
trum moves into the probed region. The fits produced 
spectral parameters a = 0.61 and /3 — 0.39. Although 
not shown here, the longitudinal data from Refs. [Sll, Isl] 
were refit using the same procedure for purposes of con- 
sistency and comparison, yielding spectral parameters 
a = 0.69 and /? = 0.30. 

Fig. |9] presents the fitted values of Tg and t; as a func- 
tion of temperature along with their respective VFT fits. 
There is excellent agreement between the two timescales 
up to the temperature T = 231 K, where the character- 
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T(K) 

FIG. 9. Characteristic relaxation time for shear (r^) and 
longitudinal (rj) degrees of freedom plotted versus temper- 
ature (left ordinate). Refitted values of the VFT parame- 
ters of the original longitudinal data are tq = 2.5 ± 1 /is, 
B = 230 ± 40 K, and To = 198 ± 4 K while the shear 
fits produced tq = 0.26 ± 0.2 /is, B = 450 ± 240 K, and 
To — 180 ± 14 K. Also shown is the decoupling parameter 
logTs — logr; as a function of temperature (right ordinate). 



istic relaxation times clearly separate. When combined 
with the observation that the HN spectral parameters 
a and /3 differ between the two degrees of freedom, we 
conclude that the shear acoustic wave dynamics differ 
significantly, at least at the lower sample temperatures 
measured, from those obtained from the earlier polarized 
ISTS data [3l|, Is^l- This conclusion, based on a larger 
data set taken with a broader range of acoustic wavevec- 
tors, supercedes the results of the previously published 
work. 

We now use our data to examine the phenomenon of a- 
scale coupling described briefly in the Introduction. Us- 
ing the fitted values of Ts and r/ as a function of temper- 
ature, we define the coupling parameter as logr^ — logr/, 
which is plotted alongside the relaxation time in Fig. |9l 
The coupling parameter is constant and essentially zero 
until it begins to decrease with decreasing temperature, 
reflecting the differing characteristic timescales of the 
shear and longitudinal relaxation as the former becomes 
slower. This provides an estimate of the MCT crossover 
temperature Tc = 231. 

We may attempt to further understand our results in 
the mode-coupling framework of Ahluwalia and Das fs^ . 
Briefly, when a calculation of a collection of hard spheres 
in a Percus-Yevick approximation is considered, the crit- 
ical length scale Lq — 27r/(7o above which propagation of 
shear waves becomes overdamped obeys the power law 
in the vicinity of the critical control parameter 



Ln = 



A 



(Ac - A)^ 



(30) 



where A — 1, 6 — 1.2, A represents the packing frac- 
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FIG. 10. Computed value of length scale limit for shear wave 
propagation as a function of temperature & = T — T' using 
the decoupling value for T' = Tc (blue circles) as the critical 
temperature and the glass transition temperature T' — Tg 
(green diamonds). For comparison, results from Ref. [30| are 
also shown as a function of packing fraction A (red squares) . 
The exponents are from the data plotted with empty symbols. 



tion, and Ac is the MCT critical packing fraction beyond 
which shear wave propagation is allowed for all length 
scales. 

Using the theoretical results of Sec. IIIIl we may at- 
tempt to deduce a lower length scale for shear wave prop- 
agation as a function of temperature by considering at 
which wavevector shear wave propagation becomes over- 
damped. For this analysis, we chose the temperature T to 
be the independent parameter, yielding a similar power 
law Lf) — A/{T — T'Y , where T' represents a critical 
temperature. Here we consider two significant tempera- 
tures for the liquid, i.e., the glass transition temperature 
Tg = 202 K [33] and the crossover temperature as 
determined from the MCT decoupling analysis above. 

Figure (TU] shows a plot of the derived upper bound, 
as deduced from Eq. [1^] as a function of the variable 
e = (T - T'). The results of Ref. ^ are also shown 



as a function of packing fraction A for comparison. A 
power law fit is shown to the four points in the vicin- 
ity of the relevant temperatures. Picking T' — Tc as 
the relevant temperature yields a fitted value of 6 = 1.2, 
which is in excellent agreement with the theoretical re- 
sult. A similar fit to the data using the literature value of 
T' = Tg produces a significantly higher value of ^ = 13.4. 
We remark that the wavelength scales Lq reached as this 
temperature is approached from above are several orders 
of magnitude larger than those corresponding to any di- 
verging structural correlation length scale. 

VI. CONCLUSIONS 



We have used depolarized impulsive stimulated Bril- 
louin scattering to measure shear acoustic waves in su- 
percooled triphenyl phosphite from 220 K to 250 K com- 
bined with previously obtained results, we are able to 
examine a frequency regime from ^10 MHz to almost 
1 GHz. Our results indicate that the shear and longitu- 
dinal spectra do not share the same spectral parameters 
a and f3. 

We also observed a-decoupling of the longitudinal and 
shear degrees of freedom, yielding an estimate of the 
mode-coupling Tc = 231 K. Using the decoupling result, 
we verified a power law for the diverging lengthscale of 
shear wave propagation as d = 1.24. A similar test using 
the literature value of Tg was not in agreement with the 
theoretical model, as should be expected since Tg is not 
a mode-coupling theory parameter. 

Further work in the study of shear relaxation in triph- 
enyl phosphite and other liquids will center on expanding 
the dynamic range of the measurements. We also note 
that a comparison with dielectric data via the model of 
DiMarzio and Bishop [1^ could also be performed if di- 
electric relaxation measurements at a similar combina- 
tion of temperatures and frequencies were carried out. 
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